TRIVIALITY OF THE QUARK PROPAGATOR IN THE 
LADDER APPROXIMATION IN QCD 



V. Gogohia 

Research Center for Nuclear Physics (RCNP), Osaka University 
Mihogaoka 10-1, Ibaraki, Osaka 567-0047, Japan 
gogohia@rcnp.osaka-u.ac.jp and gogohia@rmki.kfki.hu 



Abstract 

The validity of the ladder approximation (LA) in QCD and QED in the 
context of the corresponding Schwinger-Dyson (SD) equations and Slavnov- 
Taylor (ST) and Ward-Takahashi (WT) identities is investigated. In con- 
trast to QED, in QCD because of color degrees of freedom the summation of 
the ladder diagrams within the Bethe-Salpeter (BS) integral equation for the 
quark-gluon vertex at zero momentum transfer on account of the correspond- 
ing ST identity does provide an addition constraint on the quark SD equation 
itself. Moreover, the solution of the constraint equation requires the full quark 
propagator should be almost trivial (free-type) one, i. e. only almost trivial 
quark propagator is allowed in the LA to QCD. This triviality results in the 
fact that the standard LA ignores the self-interaction between gluons caused 
by color charges (non-Abelian character of QCD). 

PACS numbers: ll.15.Tk, 12.38.Lg 
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I. INTRODUCTION 



It is well known that the full dynamical information of any quantum field gauge 
theory such as QCD is contained in the corresponding quantum equations of motion, 
the so-called Schwinger-Dyson (SD) equations for propagators (lower Green's functions) 
and vertices (higher Green's functions) [1]. The Bethe-Salpeter (BS) type integral equa- 
tions for higher Green's functions and bound-state amplitudes [2] should be also included 
into this system. It should be complemented by the corresponding Slavnov- Taylor (ST) 
identities [1,3,4] which, in general, relate lower and higher Green's functions to each 
other. These identities are consequences of the exact gauge invariance and therefore 
"are exact constraints on any solution to QCD" [1]. Precisely this system of equations 
can serve as an adequate and effective tool for the non-perturbative approach to QCD [1,5]. 
It is important to understand however, that the above mentioned system is an infinite chain 
of strongly coupled highly nonlinear integral equations, so there is no hope for an exact 
solutions. For this reason, some truncation scheme is always needed in order to make these 
equations tractable for getting physical information from them. 

Because of its relative simplicity the most popular is the ladder approximation (LA) in 
various forms. In general, it consists of approximating the full vertices by their free pertur- 
bative (point-like) counterparts in the corresponding kernels of the above mentioned integral 
equations. Approximating in addition the full gluon propagator by its free perturbative ex- 
pression, for example in the BS-type integral equation for the vertex, one gets the so-called 
quenched LA (QLA) scheme which for the quark SD equation became known as "rainbow" 
approximation. However, very soon it was realized that QLA was too crude for QCD and 
should be improved in order to incorporate the QCD renormalization group results for the 
running coupling constant (asymptotic freedom) into the system of SD equations. The so- 
called improved LA (ILA) was proposed [6] (special dependence of the full gluon form factor 
on its momentum). It makes it possible to take into account self-interacting gluon modes 
(non-Abelian character of QCD) at the level of the full gluon propagator only, so vertices 
remain intact, i. e. point-like ones. 

The main problem of the LA is, of course, its self-consistency. For a long time, it was 
widely believed that the quark sector could be effectively decoupled from other sectors in 
QCD in the LA. For example, it is usually assumed that ghost degrees of freedom are not 
important for the above mentioned SD system of equations in this case. At first sight this is 
so indeed, since neither the quark SD equation nor the equation for the quark-gluon vertex 
depends explicitly on ghosts in the LA. However, in QCD as it was mentioned above, these 
two equations should be complemented by the corresponding ST identity which explicitly 
depends on ghost degrees of freedom (even in the LA). Unlike QED, QCD is much more 
complicated gauge theory. It contains many different sectors. Making some truncation in 
one sector, it is necessary to be sure that nothing is going wrong in other sectors since the 
above mentioned SD system of equations may remain strongly coupled even after truncation. 
The key elements relating different sectors in QCD are the above mentioned ST identities, 
playing thus an important role in the investigation of the problem of self-consistency of 
any truncation scheme. In other words, any solution to QCD, in particular LA should be 
compatible with these identities and not vice versa as it has been emphasized above. 

Especially important the problem of self-consistency of the LA becomes when it is ap- 
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plied to such essentially non-perturbative problems as quark confinement and dynamical 
chiral symmetry breaking (DCSB). Its inconsistency for the above mentioned problems was 
pointed out by Adler [7] who explicitly showed (within precisely the dynamical equations 
approach) that a more sophisticated approximation than the LA is needed to understand 
quark confinement and DCSB in QCD. The importance of the self-consistent treatment of 
DCSB within the SD system of equations in the LA was, apparently for the first time, em- 
phasized in Ref. [8]. The next interesting step in this direction was done quite recently in 
Ref. [9] where a systematic method was developed for obtaining consistent approximations 
to the SD and BS equations which maintained the external gauge invariance (by using the 
Ward-Takahashi (WT) identities for the color- singlet, but flavor-nonsiglet vertices). In con- 
trast to the above mentioned paper [9], the main purpose of this paper is to investigate the 
self-consistency of the LA by using the ST identity for the color- nonsiglet, but flavor- singlet 
quark-gluon vertex in QCD (the internal gauge invariance). Precisely this vertex is one of 
the most important key objects in QCD itself. Everybody knows that the LA is bad (for ex- 
ample, it is explicitly gauge-dependent truncation scheme), nevertheless eveybody continues 
to use it. In some sense it was justified (see, for example Ref. [10]) since up until now there 
was not an exact criterion to prove or disprove the LA in gauge theories. Here we precisely 
propose how to formulate this criterion. 

In section 2 we remind some well-known things from QED in order to better explain 
our method. In section 3 we derive the constraint equation and in sections 4 and 5 it 
is solved without and with ghosts, respectively requiring the full qurk propagator to be 
almost free one. In section 6 the flavor-nonsinglet but color-singlet axial-vector vertex is 
investigated. Some general remarks on the non-perturbative renormalization in the Landau 
gauge is presented in section 7. Our conclusions are given in section 8. 



II. QED 



Let us start with some well-known things from QED and consider first unrenormalized 
(for simplicity) SD equation in the LA for the quark (electron) propagator in momentum 
space 

r H n l 

S-\ P ) = S \p) -g 2 j j^laS(lh,D a ,(q), (2.1) 

where Sq 1 (p) = —i(p — m ) with m being the current ("bare") mass of a single quark 
(electron) and D a p(q) is the full gluon (photon) propagator in the arbitrary covariant gauges, 



D a p{q) = 



9af3 



QaQ/3 



(f 



T 



(2.2) 



and £ is the gauge-fixing parameter (£ = 0, Landau gauge). Here q = p — I is the transfer 
momentum and we assign the factor — ig^y to the point-like ("bare") veritices with the 
corresponding Dirac indices. Differentiating Eq. (2.1) with respect to p^, one obtains 

f d n l 

d^S~\p) = -r^ - d^g 2 J -^—^l a S{l)-ipD a p(q) 

= + 9 2 J j^7a[d»S(l)]^D a p(q), (2.3) 
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up to unimportant total derivative (which is assumed as usual to vanish at the ends of 
integration). This is the differential form of the quark SD equation (2.1) relevant for further 
discussion in the next section as well. 

Let us consider now the Bethe-Salpeter (BS) integral equation in the LA for the quark- 
photon vertex T^(p, q) (q — p' — p) at zero momentum transfer q — 0, 

r M (p,0) = -il, -9 2 J -^ la S(l)T,(l,0)S(l)^D a ^q). (2.4) 

The corresponding WT identity in QED, q^T^p + q,p) = S" 1 ^ + q) — S~ l {p), provides an 
exact solution for T M (p, 0), namely 

T,(p,0)=d,S-\p). (2.5) 

Substituting it into the Eq.(2.4) and using the obvious identity 9 M 5' _1 (p) = 
— S~ 1 (p)[d fl S(p)]S~ 1 (p), one finally obtains 

r d n l 

d,S-\p) = ^ + g 2 j j^l a [d,S(l)h,D a ,(q). (2.6) 

This is the quark SD equation obtained from the BS integral equation for the corresponding 
vertex on account of the WT identity (2.5). Comparing it with that of Eq. (2.3) (second 
line), one immediately concludes that they are identical. The renormalized version of Eq. 
(2.3) is obtained by the multiplication of its inhomogenious term by the quark wave function 
renormalization constant Z 2 , while the renormalized version of Eq. (2.6) is obtained by 
the multiplication of its inhomogenious term by the vertex renormalization constant Z\. 
Because of the WT identity (2.5) these constants are equal to each other [Z\ = Z 2 ) and 
again renormalized versions of these equations coincide in QED. So the BS equation at 
zero momentum transfer and on account of the corresponding WT identity (which is a 
consequence of the exact gauge invariance) does not provide an additional constraint on the 
quark SD equation in QED. In other words, the summation of the ladder diagrams within 
the corresponding BS integral equation at zero momentum transfer on account of the WT 
identity is consistent with the electron (quark) SD equation itself. 



III. QCD 

The situation in QCD is completely different and is much more complicated because of 
color and ghost degrees of freedom which precisely provide an addition constraint on the 
solution to the quark SD equation in the LA by using the same method as in previous section 
for QED. It is convenient to begin our analysis from the quark SD equation in QCD which is 
just the same as in QED (2.1) with only trivial replacement g 2 — > g 2 F because of color group 
factors. Here g F = g 2 Cp and Cf is the eigenvalue of the quadratic Casimir operator in the 
fundamental representation (for SU(N), in general, C F = {N 2 - 1)/2N = 4/3, iV = 3). Thus 
the renormalized version of the quark SD equation (2.1) is as follows 

S- 1 (p) = Z 2 S 1 (p)+zt(p) 

= Z 2 S^(p) -g 2 F J j--l a S{l)^D aP {q), (3.1) 
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where the renormalized quark propagator and the full gluon propagator (2.2) are related to 
the unrenormalized ones as 



S(p) = Z 2 S(p), 
D a p(q) = Z v D a/3 (q) 



(3.2) 



with the renormalization of the coupling constant is beeing determined by the corresponding 
combination of the above and below introduced renormalization constants, Z's [4]. The 
differential form of the renormalized quark SD equation becomes (compare with Eq. (2.3)) 



d l ,S-\p) = -Z 2 i 1 ^ + d^{p) 



-Ziiip + g 2 F 



d n l 



■la[d,,S{l)YlpD a p{q). 



(3.3) 



Let us introduce now the renormalized ST identity [1,3] 



7-1 



b(k 2 ; Oj = yZg T a - B a (p, k; £)J S~\p + k) - S'^p) [Z^T a - B a ( P} k; £) 

(3.4) 



where the corresponding renormalized quark-gluon vertex, the ghost self-energy and the 
ghost-quark scattering kernel[] are related to their unrenormalized counterparts as follows: 



r„(p, k) = Z 1 1 f lx {p, k), 

b(k 2 -0 = Z 9 b(k 2 ;0, 
B(p,k;0 = Z B B(p,k;0. 



(3.5) 



For furture aim we have introduced an explicit dependence on a gauge fixing parameter £ 
into the ghost degrees of freedom. Let us note that in the ST identity (3.3) the transfer 
momentum runs through the ghost degrees of freedom. The following important relation 
between renormalization constants holds [1] 



Z\Zb — Z2Z q . 



(3.6) 



It is worthwhile noting that for the sake of furture convenience instead of the ghost field 
renormalization constant and the "renormalization" constant Y for the amplitude H(p, k) = 
1 — B(p, k), defined in Ref. [1], we have introduced the ghost self-energy remormalization 
constant Z g and the "renormalization" constant Zb for the ghost-quark scattering kernel 
itself. 



1 For its skeleton expansion see, for example papers [11,12] and references therein. In the LA only 
the first term of this expansion survives with the replacement of the corresponding full ghost- and 
quark-gluon vertices by their free perturbative (point-like) counterparts. However, there is no need 
to exploit the LA explicit expressions for quantities defined in (3.8) and (3.9) (see below) within 
our approach. 
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Differentiation of the renormalized ST identity (3.3) with respect to and then setting 
k = 0, yields^ 

f -(p, o) [z; 1 + 5(o; o] = z^T^s-'ip) - n»(p, o; o - *>; fls- 1 ^) + s- x (p)*>; 0, 

(3.7) 

where we introduced the following notations 



and 



dk^ 



(3.8) 

fc=0 



n«(p,0;0 = S B (p,0;0^" 1 (p)- ( 3 - 9 ) 

Here a few remarks are in order. Though the regular dependence of the ghost-quark scat- 
tering kernel B a (p,k;^) on the ghost self-energy momentum k is preserved by the general 
Taylor's result [14] 

B a (p,0;0 = £F a (p) = at ^ = 0, (3.10) 

which shows that this kernel exists at small k in any covariant gauge and vanishes only in the 
Landau gauge £ = 0, nevertheless the regular dependence of the ghost self-energy itself on 
its momentum is not so obvious. However, here we treat the ghost-self-energy as a regular 
function of its momentum since a singular dependence (which, in principle should not be 
excluded a priori) requires completely different investigation and is left for consideration 
elsewhere. The constraint equation in this case also will be completely different from that 
obtained below. 

The renormalized version of the BS-type integral equation for the color non-singlet vertex 
at zero momentum transfer is 

r>, 0) = -ZtvypT* -~g 2 j —-- la T b S{l)T^ 0)S(l)^T b D a p( q ). (3.11) 

Substituting (3.7) into the BS-type integral equation (3.11), one obtains 
Z^T^S-^p) - ft«(p, 0; - 0S-\p) + S~ l (p)^(p; = -Z X \Z? + 6(0; £)]i 7 ^° 

+ ?Z B ^T b T a T b J -^l a [d,S{l)] lp D aP {q) 



r d n l ~ 
f J j-~l a T b S{l)W^ 0; 0S(lh,T b D a ,(q) 



(2tt 
d n l 



/d ~ 
7^^T b ^(l-,0S(lh P T b D aP (q) 



(2vr 
d n l 



/(I ~ 
j^l a T b S(l)^(l; Z)jpT b D aP ( q ), (3.12) 



2 It is assumed that possible unphysical kinematic singularities have been already removed from 
the vertex by means of Ball and Chiu procedure [13]. 
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where the obvious identity, d ll S~ 1 {p) = — S~ 1 (p)[d fl S(p)]S'' 1 (p), has been already used. Let 
us now use the commutation relation between color matricies [T a , T b ] = if a bcT c , where f a b c 
are the antisymmetric SU(3) structure constants with non-zero values, given for example in 
Ref. [15]. Then one obtains, 

T b T a T b = p p _ }_c A ]T a , (3.13) 
2 

where C F is the above mentioned eigenvalue of the quadratic Casimir operator in the fun- 
damental representation while Ca is the same but in the adjoint representation, C\ = N for 
SU(N) (N = 3 for QCD). So from this relation and on account of the renormalized version 
of the differential form of the quark SD equation (3.3) (because of g 2 Cp = g F ) and using 
finally an important relation (3.6), one arrives at 

- H>, 0; - *J(p; 0S~\p) + S-\p)%(p; = -Z^O; 

1 r d n l 

- ^C A T a fZ^j ^7«[W)]7/J^(?) 

r d n l ~ 
+ ~ 9 J j2^ laTbSma ^ ° ; 0S(lh,T b D a ,(q) 

r d n ] ~ 
-f\ j^l«T b ^(l-,OS(lhpT b D a p( q ) 

r d n l ~ 

+ 9 2 J J^laT b S{l)^{l-i) lp T b D aP {q). (3.14) 

This is a renormalized version of the general constraint equation which relates ghost degrees 
of freedom to those of quark ones in a new manner within the LA. In principle, the solution 
(if any) of the quark SD equation (3.1) itself should be compatible with this constarint 
equation. At first sight it is too complicated integral equation, but nevertheless its analysis 
is rather simple. Here let us note that constraint (3.14) is slightly simplified in the Landau 
gauge £ = since the composition (3.9) vanishes due to the above mentioned Taylor's general 
result (3.10). 



IV. SOLUTION WITHOUT GHOSTS 

As it has been already mentioned in the Introduction, it is widely believed that the 
ghost degrees of freedom are not important for the LA to covariant gauge QCD. If this is 
so indeed, let us omit them from the general constraint equation (3.14) by "hand", which is 
equivalent to consider QED-type ST identity (3.4) (apart from the color group generators) 
from the very beginning. Apart from neglecting corrections to the vertex, this is the standard 
procedure in the LA to covariant gauge QCD. It has been already done in order to proceed 
from the SD equation for the vertex to its BS-type counterpart (3.11) which allows precisely 
to explicitly sum up ladder contributions (see, for example Ref. [5]). Doing so in (3.14) and 
on account of Eq. (3.3), one obtains 

r d n l ~ 

J 7^^[d,S(l)]^D a ^q) = ~g F %iE(p) = 0, (4.1) 
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which implies general solution is to be i£(p) = im c , where m c is the constant of integration 
of the dimension of mass. From the quark SD equation (3.1) (first line) 

S- 1 (p) = Z 2 S 1 (p)+zt(p) (4.2) 

with Sq 1 ^) = —i(p — m ) and on account of the above mentioned trivial solution to the 
quark self-energy, it finally follows 

Sip) = J^, (4.3) 
p — m 

where Z 2 l = Z 2 and rh = m + Z 2 m c . Thus the solution of the constraint equation (3.14) 
with ghosts omitted by "hand" in the LA to covariant gauge QCD requires that the quark 
propagator should be almost free one apart from the redefinition of the quark mass. In 
other words, there should be no nontrivial solution to the quark SD equation in the LA in 
this case. So the main question arises, namely may ghosts cure this problem or not? The 
final answer is, of course, not (see next section). If they would cure this problem it would 
mean that the quark propagator (via the general constraint equation (3.14)) would explicitly 
depend on ghost degrees of freedom. However, this is impossible in QCD [1] where ghosts 
contribute only into the closed loops, so nothing can explicitly depend on them exept of the 
contributions to the renormalization constants, constants of integration, etc. Precisely the 
above mentioned problem will be investigated below in more sophisticated fashion. 



V. SOLUTION WITH GHOSTS 

In order to investigate the general constraint equation (3.14) in more sophisticated way, 
let us write down formal BS-type integral equations for the compositions defined in (3.8) and 
(3.9). Like the quark-gluon vertex, these objects are vector quantities at zero momentum 
transfer. They satisfy formally absolutely the same BS-type integral equations in the LA 
as the above mentioned quark-gluon vertex at zero momentum transfer (3.11). Thus they 
are the LA equations for (3.8) and (3.9) in the same way as (3.11) is the LA to the exact 
BS integral equation for the quark-gluon vertex at zero momentum transfer. The one of 
the differences is that there are no point-like counterparts of these quantites in QCD. But 
the main difference is that, in contrast to the BS-type integral equation for the quark-gluon 
vertex (3.11), these equations are completely auxiliarly. They have no independent role, 
the main purpose to use (or to postulate) them is to show (as was mentioned above) that 
nothing explicitly depends on ghost degrees of freedom in QCD. So in complete analogy 
with (3.11), the formal homogenious BS-type equation in the LA for the composition (3.9) 
looks like 

r d n l ~ 

n>, 0; = -f J j^7«T b S(l)n;(l, 0; i)S{l) lp T b D aP {q). (5.1) 

For the composition defined in (3.8) it is convenient to write down the formal BS-type 
integral equations for its left and right hand side combinations, ^^(p;^) = ^(p; C)S~ l {p) 
and #° (r) (p;£) = 5 _1 (p)^(p;0, respectively 
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r d n l 

<°(P;0 = ~f\ ^7 Q T^(/)^f)(/;O^(07^ J D a/3 (g) 
r c\ n ] ~ 

*? r) (p;0 = -£ 2 / ^7«T*5(0^(i; 05(07^^(9). (5.2) 

Substituting (5.1) and (5.2) into the constraint equation (3.14), we are left with 

1 r d n l 

- Z 1 b(0; Oh, = ~ 2 C A fZ B ' J —^S(l)hf3D af3 (q). (5.3) 

From Eq. (3.3), however it follows 

r d n l ~ 

J (2^^[W)]7/j£>au(?) = ^ 2 ^iS(p), (5.4) 

so Eq. (5.3) by denoting &i(0;f) = 2Z 1 Z B (C F /C A )b(0;£) = 2Z 2 Z g (C F /C A )b(0;^, because 

of (3.6), becomes — fy M 6i(0; £) = <9 M i£(p). Its general solution is, = — ip&i(0;£) + im c , 

where m c is the constant of integration of the dimension of mass. From the quark SD 
equation (4.2) it finally follows 

S(p) = i (5.5) 
p — m 

where 

Z 2 l = Z 2 (l +2Z 9 (C F /C A )b(0;0), 

m = Z 2 (m c + Z 2 m ). (5.6) 

Thus the use of the auxiliarly BS-type integral equations (5.1) and (5.2) allows one to 
treat ghost degrees of freedom in more sophisticated way than simply omit them by "hand" 
(as in previous section), but, nevertheless the quark propagator in the LA again remains 
trivial one (5.5) apart from the redefinitions of the quark mass and the quark wave function 
renormalization constant with the help of ghost self-energy at zero point (compare with 
(4.3)). Nothing explicitly depends on ghosts as it should be indeed in QCD. It is well known 
that in the Landau gauge £ = the quark wave function is not renormalized, i. e. Z 2 — 1. 
Whether the redefined renormalization constant Z 2 is also will be equal to unity in this gauge 
depends on the solution of the SD equation for the ghost self-energy. This investigation as 
well as the solution of the quark SD equation itself in order to determine a possible value of 
the constant of integration m c is, obviously, beyond the scope of this paper and is left for 
consideration elsewhere (see also section 7). 



VI. QCD. FLAVOR-NONSINGLET AXIAL- VECTOR VERTEX 

It is convenient to continue our discussion of the situation in QCD with the ST identity 
for the flavor-nonsinglet (but color- singlet), axial- vector vertex. It is not complicated by 
the unknown ghost contributions. In QED there is no sense to consider axial-vector vertex 
because there exists only flavor-singlet channel which, obviously, will be always complicated 
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by the axial anomaly [16]. Contrast to QED, in QCD there exists flavor-nonsinglet channel 
which is free from axial anomaly, but a new interesting problem occurs with this identity. 
Also not loosing generality, let us consider the chiral limit (m = 0) of this identity since 
the nonchiral case, evidently, can not change conclusions drawn below. 

Let us consider the flavour non-singlet, axial- vector WT identity in the chiral limit 

<^n> + ?,p)= [jj{^S- 1 (p) + S-\p + q) l5 }, (6.1) 

where q = p' — p is the momentum transfer (the external momentum) and the quark prop- 
agator is given by 

-tS(p)=pA(p 2 ) + B(p 2 ), (6.2) 

so its inverse is expressed as 

{-iS{p)y l =pA{p 2 )-B{p% (6.3) 

with 

A{p 2 ) = A(p 2 )D~\p 2 ), 
B(p 2 ) = B(p 2 )D-\p 2 ), 

D{p 2 )=p 2 A 2 {p 2 )-B 2 {p 2 ). (6.4) 

In connection with (6.1), one has to point out that, in general, X 1 is a SU (Nf) flavour matrix 
and, in the massless case, the quark propagator is proportional to the unit matrix in the 
flavour space. 

It is well known that if dynamical chiral symmetry breaking (DCSB) at the quark level, 
implemented as 

{S- 1 (p), l5 } + = t l5 2B(p 2 )^0 (6.5) 

takes place, then the axial-vector vertex should have a pole corresponding to a Goldstone 
state. Indeed, from (6.5) it follows that T\ {jp + q,p) has a pseudoscalar pole at q 2 = 
(dynamical singularity which determines Goldstone state) if and only if the double of the 
dynamically generated quark mass function 2B(p 2 ) in (6.5) is nonzero and vice versa. In 
order to self-consistently untangle regular and pole parts in (6.1), it is necessary to perform a 
procedure described in detail in our paper [12]. It is a similar to that of the above mentioned 
Ball and Chiu procedure [13] to remove kinematic (unphysical) singularities from the vertex. 
As a result of this, dependence on the two arbitrary form factors occurs in the regular part 
of this vertex at zero momentum transfer q — 0, namely 

rt>,p) = (jj [-id,S-\p) + A»] 75 , (6.6) 

where 

A„(p) = iMp 2 ) - PIMP 2 )- (6.7) 
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The corresponding BS integral equation becomes equation for these form factors and there- 
fore no statement about justification of the LA can be deduced. Indeed, let us show this 
explicitly. 

The BS integral equation for the regular part of the axial- vector vertex at zero momentum 
transfer is (here and below in this section again q—p — l) 

yj 757, -glj j^l«Sm% l)S(lh,D afS (q), (6.8) 

where summation over color indices is already done, i. e. g 2 F = g 2 Cp (see Eq. (3.1)). 

Substituting relations (6.6-6.7) into this equation and on account of the above used 
identity, one obtains 

[d^S^ip) + iA M (p)]7 5 = «757 M 

g2p I (^^[ 9 ^( / )]^' 1 (0755(/)7/ 3 ^/ 3 (g) 
r d n l 

" ^ J J^^S^A^S^pD^q). (6.9) 
Let us now introduce the following notation 

S-\ph 5 S(p) = -75 + A(p) 75 . (6.10) 

where, obviously A(p) can be given in terms of the quark propagator form factors (6.4). 
Multiplication of Eq. (6.9) by 75 from right and on account of (6.10), yields 



+ 



d fl S- 1 (p)+tA fl (p) = -t llA 

d n l 



■7a[<9 M 5(/)]7 /3 D a ^(g) 



(27T) 

r d n l 

~ 91 J 7^7a[^(0]A(07^(?) 
r d n l 

~ igl j j^l a S(l)A,{l) lb S{l) lpl ,D a ^q). (6.11) 

Using futher the QCD quark SD equation (3.3), one finally obtains 

<A M (p) = -9 2 F j j^l a [d,S{l)]A{l) 1(3 D aP {q) 

-SfJ ^laS(l)A,(l) l5 S(lh^ 5 D a(S (q), (6.12) 

where, recalling A M (p) is given in (6.7). This is the inhomogenious BS integral equation 
for unknown form factors Rq and Ru, so, as was mentioned above, no statement about 
justification of the LA can be deduced. Evidently, this is a result of the fact that the exact 
expression for the longitudinal part of color-singlet but flavor-nonsinglet axial-vector ST 
identity is known up to two arbitrary form factor even at zero momentum transfer. Also 
obviously, that one can derive the similar BS equation for these form factors in some other 
approximation for the BS scattering kernel and again no statement can be achived in favour 
of this or another truncation scheme. 
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VII. REMARKS ON THE NON-PERTURBATIVE RENORMALIZATION IN 

THE LANDAU GAUGE 



In this section we are not restricted to the LA and are going to make a few general remarks 
concerning the non-perturbative renormalization in the Landau gauge. It is well-known that 
in QCD the renormalization program at zero point is greatly simplified in the Landau gauge. 
The ghost degrees of freedom, such as the ghost-gluon vertex and the ghost-quark scattering 
kernel, do not renormalize at zero point in the Landau gauge. Their renormalization con- 
stants (defined at zero point) become trivial unites in this gauge. Summary of these results 
one can find in Ref. [1]. This fails if the renormalization is performed at nonzero euclidean 
mass scale parameter ((— /z 2 ), renormalization point) because of the possible infrared (IR) 
singularities due to exchange of massless gluons. However, as was mentioned above, the IR 
singular behaviour of the ghost self-energy requires special treatment and will be considered 
elsewhere .0 

Let us consider now the exact ghost propagator and its self-energy which are related to 
each other as follows 

where, for simplicity's sake, we suppressed the color group indicies. Introducing the renor- 
malized ghost propagator 

G(k) = Z G G(k), (7.2) 
from (7.1-7.2) and (3.5), one obtains 

Zg = Z- 1 (7.3) 

and 

e <*> = w^moY {7A) 

where we define as usual (see, for example, Ref. [4]) 

6 B (* a ;O = 6(* 2 ;O-S(0;O (7-5) 

and 

Z- 1 + 6(0;0 = l- (7.6) 



3 The treatment of ghost degrees of freedom in the Landau gauge in Ref. [17] is not completely 
self-consistent. It was not investigated whether the enhancement of the ghost propagator in the 
IR was compatible with Taylor's general result (3.10) for the ghost-quark scattering kernel which 
was simply omitted by hand in their truncation scheme. 
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In a similar way let us define the regularized quark-ghost scattering amplitude (3.5) as 
follows 

B R (p,k;t) = B(p,k;t)-B(0,0;t) (7.7) 

and 

Z B l - £(0,0; = 1, (7.8) 

so that 

Z B 1 -B(p,k;0 = l-B R ( P ,k;0 (7.9) 

From Taylor's general result (3.10) and on account of (7.8) then it follows that in the Landau 
gauge £ = 0, one has 

Z B = 1, (7.10) 

and consequently from (3.6) it follows 

Z 1 = Z 2 Z g . (7.11) 

It is worthwhile noting that the renormalization constant Zb coincides with the renormal- 
ization constant Y defined in Ref. [1] for the amplitude H(p,q) = 1 — B(p,q), i.e. Zb = Y 
in the Landau gauge £ = at least. 

Let us now briefly discuss a possible interesting feature of the renormalization at zero 
point in the Landau gauge. Indeed, from (7.11) it follows that if the ghost self-energy 
renormalization constant Z g (defined at zero point) would may equal to unity as Zb (7.10) 
then the QED-type relation between the QCD renormalization constants (defined at zero in 
the Landau gauge) will be recovered, namely (7.11) would become 

Z x = Z 2 . (7.12) 

But from (7.6) then it follows that the ghost self-energy at zero point in the Landau gauge 
should become zero, i.e. 6(0; £ = 0) = 0. This would be possible only within the nonpertur- 
bative treatment of the corresponding ghost self-energy SD equation. The problem is that 
the perturbative calculation of this renormalization constant within the LA to the above 
mentioned SD equation suffers from the IR singularities due to exchange of massless glu- 
ons. In this case this would mean that the ghost propagator in the Landau gauge does not 
renormalize at zero point as well as the ghost-gluon vertex and ghost-quark scattering kernel 
[1]. This would be an example of the nonperturbative finite calculation of the corresponding 
renormalization constant and deserves furture investigation elsewhere. 

A few concluding remarks are in order. The SD ghost self-energy equation is complicated 
nonlinear integral equation relating three various Green's functions to each other. There are 
no doubts that this equation may have regular as well as singular solutions even in the LA. 
However, it is necessary to clearly understand that the disappearance of the ghost-self energy 
at zero point in the Landau gauge is completely nonperturbative effect. One can not obtain 
this result calculating only finite number of perturbative corrections. Only summing up 
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infinite number of perturbative diagrams within the corresponding SD equation (i.e. going 
beyond the perturbation theory), one may hope to achieve this goal. Moreover, it depends 
heavily on the nature of the approximation scheme, i.e. not any approximation scheme 
may satisfy this exact result. In other words, the above mentioned possible zero ghost self- 
energy in the Landau gauge is an example of the nonperturbative finite calculation of the 
corresponding renormalization constants at zero point since the perturbative renormalization 
at this point suffers from the IR singularities due to exhange of massless gluons. Precisely 
for this reason, the perturbative renormalization is always performed at nonzero space-time 
momentum k 2 = —fi 2 in order to avoid (not to solve, of course) these difficulties. At the 
same time, as it seems to us, the nonperturbative renormalization will make it possible 
to perform finite calculations of the corresponding renormalization constants at zero point 
without any problems, i. e. only beyond the perturbation theory one is able to renormalize 
theory at zero point at least in the Landau gauge. This emphasizes once more the special 
role of the Landau gauge in QCD. 

VIII. CONCLUSIONS 

In summary, we have formulated the method how to prove or disprove the LA in gauge 
theories such as QED and QCD. In QED the summation of the ladder diagrams within the 
BS integral equation for the quark-photon vertex at zero momentum transfer on account of 
the corresponding WT identity does not provide an addition constraint on the solution to 
the quark SD equation itself. In other words, there is no criterion to prove or disprove the 
use of the LA in QED. In contrast to this, in QCD because of color degrees of freedom the 
summation of the ladder diagrams within the BS integral equation for the quark-gluon vertex 
at zero momentum transfer on account of the corresponding ST identity does provide an 
addition constraint on the solution to the quark SD equation itself. Moreover, the solution 
of the constraint equation (3.14) requires that the full quark propagator should be almost 
trivial (free-type) one (4.3) or (5.5), i. e. there is no nontrivial quark propagator in QCD 
in the LA. In other words, there is no running quark mass in the LA to QCD as well. This 
is our main result. It does not depend on how one treats ghost degrees of freedom in the LA 
to covariant gauge QCD, omitting them by "hand" in the general constraint equation (3.14) 
or in more sophisticated fashion by using auxiliarly BS-type integral equations (5.1) and 

(5.2) there. Let us underline once more, that the existence of the constraint equation (3.14) 
in the LA to QCD in any gauge is only due to color degrees of freedom and not ghosts. Let 
us make one thing perfectly clear. The quark SD equation, remainning a nonlinear integral 
equation even in the LA, may have formally a number of non-trivial, approximate solutions. 
However, non of them (analytical or numerical) will satisfy the constraint which comes from 
the exact ST identity, so all results based on these solutions should be reconsidered. 

In noncovariant gauges [18] ghosts do not contribute from the very beginning, i. e. the 
ST identity (3.4) becomes of the QED-type apart from the color degrees of freedom. This 
means that the general constraint equation (3.14) automatically becomes (4.1). Thus from 
our consideration, one easily can conclude that the quark propagator remains trivial one 

(4.3) in the LA to noncovariant gauge QCD as well (let us remind that we neither use the 
explicit expression for the full gluon propagator (2.2) nor specify the full gluon form factor 
there). The physical reason of the discovered triviality of the quark propagator in the LA 
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is, of course, that the LA in the quark sector ignores the self-interactions of gluons caused 
by color charges (non-abelian character of QCD). On the other hand, there is no doubt that 
the LA (especially in the ILA form) is legitimated to use for the ultraviolet (UV) region in 
QCD because of asymptotic freedom. The problem how to correctly go beyond the LA in 
order to describethe IR region in QCD becomes inevitably important. It is perfectly clear 
now that we can not use the LA (in QLA or even in ILA forms (since as it was mentioned 
above we did not specify the full gluon form factor)) in the entire energy- momentum range 
in QCD. The ILA itself is not sufficient to approximate the IR region in QCD. The nontrivial 
generalization of the point-like vertices is needed for this purpose. We have investigated the 
self-consistency of the LA in QCD. It is almost obvious that any truncation scheme, for 
example such as planar, 1/N C limit, etc maight and should be investigated in the same way. 
As it is emphasized in our paper, this is important in gauge theories. 

We didnot discuss the flavor-nonsinglet (but color-singlet) vector vertex to which external 
gauge fields such as W and Z bosons could couple. Evidently, the situation for this vertex 
is completely analogous to that of QED. In other words, contrast to the axial-vector vertex 
considered in Section 6, here the use of the LA is justified as in QED. First it was explicitly 
shown in Ref. [9] though in the different way. However, in the light of our result a general 
question is in order. How is it possible to use nontrivial solutions to the quark SD equation 
in the LA in the QCD boun-state problems when one knows that it becomes trivial (almost 
free) in the scattering problems? At long last, the quark SD equation is the same for both 
problems. At least in one case, namely in the constituent quark model (CQM) the use 
of the LA is justified for both problems. Indeed, from a QCD theoretical field point of 
view, apart from other simplifications, the CQM mainly is nothing but an approximation 
of the full quark by the constituent quark propagator (which coincedes with (5.5) due to 
replacement m — > m g , where m q is the constituent quark mass) and the full constituent 
quark-gluon vertex by point-like one. From our results then it follows that in CQM the 
quark propagator and quark-gluon vertex is in one-to-one correspondence with each other. 
Precisely, this self-consistent relation between them explains the success of this calculation 
scheme despite its simplistic structure. 

In conclusion, let us make a few technical remarks. In QED expansion in powers of the 
external momentum makes no great sense since there is no stable bound-states (positronium 
is unstable). In contrast to this, in QCD the same expansion makes precisely great sense 
because of the existence of the Goldstone sector there. Many important physical quantities 
such as scattering lenght, pion charge radius, etc are defined as coefficients of the pion form 
factor expansion in powers of external momentum. The chiral perturbation theory [19] and 
its counterpart at the fundamental quark level [20] are also based on these expansions. Thus 
in QCD zero external momentum transfer has a physical meaning as relating directly to 
various physical observables. 

The author would like to thank S.Adler for correspondence and A.T.Filippov and A. A. 
Slavnov for interesting remarks. He is also grateful to F.Schoberl and W.Lucha for critical 
discussions and useful remarks on first stage of this investigation. 
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